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Abstract 



We compute the probabilty for the processes A — > qq,gg via vacuum 
polarization in the presence of a classical space-time dependent non-abelian 
field A by applying the background field method of QCD. The processes we 
consider are leading order in gA and are simillar to ^ — > e'^e" in QED. Gluons 
are treated like a matter field and gauge transformations of the quantum 
gluonic field are different from those of the classical chromofield. To obtain 
the correct physical polarization of gluons we deduct the corresponding ghost 
contributions. We find that the expression for the probability of the leading 
order process A ^ gg is transverse with respect to the momentum of the 
field. We observe that the contributions from higher order processes to gluon 
pair production need to be added to the this leading order process. Our 
result presented here is a part of the expression for the total probability for 
gluon pair production from a space-time dependent chromfield. The result 
for qq pair production is similar to that of the e~^e~ pair production in QED. 
Quark and gluon production from a space-time dependent chromofield will 
play an important role in the study of the production and equilibration of the 
quark-gluon plasma in ultra relativistic heavy-ion collisions. 

PACS: 12.38.Aw; 11.55.-q; ll.15.Kc; 12.38.Bx 
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I. INTRODUCTION 



Over the years there have been several investigations on the production of charged par- 
ticles from a classical electro-magnetic field, a phenomenon which was discovered nearly 
fifty years ago by Schwinger By now, the production of electron-positron pairs from 
the abelian field is extensively studied both theoretically and experimentally [|2[. The sub- 
ject of quark/anti-quark and gluon pair production from the non-abelian field is relatively 
new and is not fully solved. It might be important for the production of the quark-gluon 
plasma (QGP) in the laboratory by high energy heavy-ion collisions. Lattice QCD predicts 
the existence of such a state of matter at high temperatures (~ 200 MeV) and densities 
01 . In high energy heavy-ion collisions at RHIC and LHC the receding nuclei might 
produce a strong chromofield which would then polarize the QCD vacuum and produce 
quark/ ant i- quark pairs and gluons. These produced quarks and gluons collide with each 
other to form a thermalized quark-gluon plasma. The space-time evolution of the quark- 
gluon plasma in the presence of a background chromofield is studied by solving relativistic 
non-abelian transport equation of quarks and gluons with all the dynamical effects taken 
into account. As color is a dynamical variable in the non-abelian theory, the relativistic 
non-abelian transport equation for quark and gluon is 

p,d^ + gQ^F^ya^ + gr^V^Ayd'o] f{x,p, Q) = C{x,p, Q) + Q). (1) 

Here f{x,p, Q) is the single-particle distribution function of the parton in the 14 dimensional 
extended phase space which includes coordinate, momentum and color in SU(3). The first 
term on the LHS of Eq. (Q) corresponds to the usual convective flow, the second term is 
the non-Abelian version of the Lorentz-force term and the third term corresponds to the 
precession of the color charge, as described by Wong's equation C and S on the RHS 
of Eq. dll) are the collision and the source terms, respectively. Note that there are separate 
transport equations for quarks, antiquarks and gluons The source term S contains 

the detailed information about the production of quarks and gluons from the chromofield. 
It is defined as the probability W for the parton production per unit time per unit volume 
of the phase space. Hence the space-time evolution of the quark-gluon plasma in a high 
energy heavy-ion collision crucially depends on how quarks and gluons are produced from 
the chromofield ||^. 

The production of fermion pairs from the field is studied in two different cases: 1) 
from a constant, uniform electric field and 2) from a space-time dependent field. The 
fermion pair production from a constant, uniform field is computed by Schwinger. This is 
an exact one- loop nonperturbative result This result can also be understood in terms 
of a semiclassical tunneling across the mass gap 0. However, a space-time dependent 
field A can directly excite the negative energy particles to levels above the mass gap, by a 
perturbative mechanism A qq, without recourse to any tunneling or barrier penetration 
mechanism The particle production from this mechanism is important in high energy 
heavy- ion collisions. As shown by numerical studies [^], the chromofield acquires a space- 
time dependence as soon as it starts producing partons and so Schwinger's non-perturbative 
formula for particle production from a constant field is not applicable. In this case the 
treatment of parton production from a space-time dependent chromofield is necessary. It 
can be mentioned here that, so far the transport equation is solved with parton productions 
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from a constant field taken into account because gluon pair production from the space- 
time dependent chromofield is not computed. Aim of this paper is to calculate the probability 
for the process A ^ gg which is similar to the process A qq. Quark and gluon production 
from a space-time dependent chromofield is needed to study the production and equilibration 
of a quark-gluon plasma in ultra relativistic heavy-ion collisions at RHIC and LHC 

The production of qq pairs from a space-time dependent non-abelian field is almost the 
same to that of the production of e~^e~ pairs from the abelian field. This is because, apart 
from the color factors, the interaction of the quantized Dirac field with the classical potential 
is the same in both cases. All the methods used to obtain the probability for the production 
of the e+e~ pair from a space-time dependent Maxwell field can be applied to obtain the 
probability for the production of the qq pair from a space-time dependent Yang-Mills field. 
On the other hand, the production of gluon pairs from a space-time dependent Yang-Mills 
field is not straight forward and there is no counter part to this in the abelian theory. In 
contrast to abelian theory the quantized Yang-Mills field interacts with the classical non- 
abelian potential. Because of this interaction there is gluon production from the QCD 
vacuum in the presence of an external chromofield. This phenomenon is absent in the 
Maxwell theory. 

We again mention here that the fermion pair production from the space-time dependent 
field is studied in the literature 0,^,0 but gluon production from space-time dependent 
chromofield is not studied so far. This is because a consistent theory involving the interaction 
of the gluons with the classical chromofield is not available in the conventional theory of 
QCD. In the case of fermions one knows exactly the theory for the interaction of fermions 
with the classical field. This is given by the Dirac equation (see section II). However, a 
consistent theory of gluons in the presence of external classical chromofield is not available 
in the conventional QCD. This problem is addressed properly within the background field 
method of QCD which was introduced by DeWitt and 't Hooft. We will briefly mention 
the main differences between the background field method of QCD and the conventional 
method of QCD in section III. Unlike in conventional QCD, the Feynman diagrams obtained 
in the background field method contain classical chromofields and gluons. We will use the 
Feynman rules derived in the background field method of QCD to obtain the probability 
for the processes A — > qq, gg via vacuum polarization. The study of the production and 
evolution of the quark-gluon plasma at RHIC and LHC by solving the relativistic non-abelian 
transport equation with parton production from a space-time dependent chromofield will 
be undertaken in the future. 

This paper is organized as follows. We compute the probability for the process A ^ qq 
in section II. The probability for the process A ^ gg is computed in section HI. A brief 
description about future research and conclusions can be found in section IV. 



II. PROBABILITY FOR THE PROCESS Aql ^ QQ 

The production of qq pairs from a non-abelian field via vacuum polarization is simillar to 
that of the production of e~^e~ pairs from the abelian field. This is because the interaction 
lagrangian of the quantized Dirac field with the classical gauge potential is similar in both 
the cases. The amplitude for the lowest order process A —>■ e+e~ which contributes to the 
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production of e~^e pair from a space-time dependent classical abelian field is given by 

in 

M =< A;i,A;2|^^^^|0 >= -ieu{ki)^f,A^'{K)v{k2). (2) 

In the above expression ki, k2 are the four momenta of the produced electron and positron 
respectively and A{K) is the Fourier transform of the space-time dependent field A'^{x) with 
K = ki + k2. 

The probability for pair production is: 

^ f f f (27r)^ - fci - k2) T (3) 

where 

r = s,p,„|Mp. (4) 



Performing the spin sum one obtains [11 



»'iii- = ?/, rf*A'(l-^) (l + '^-^)mK)\'-\B{K)\\ (7) 



- f t>.„, (1 - ||)"V + H) ■ MK)K . A-) 

• (5) 
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In the above expression a = |^ is the coupling constant. Using 

F^'^{K) = -iiK^A" - K'^A^') (6) 

one obtains 

This is the probability for the production of the e'^e~ pair from a space-time dependent 
abelian field. This result was obtained for the first time by Schwinger (see Eq. (6.33) in |jl|). 

Now we proceed to compute the probability for the production of a qq pair from a space- 
time dependent non-abelian field. The amplitude for the lowest order process A —>■ qq which 
contributes to qq pair production from a space-time dependent classical non-abelian field 

M = igu\k,)'j^^A''^{Ky{k2). (8) 

In the above expression ki and ^2 are the four momenta of the quark and anti-quark respec- 
tively and A"'^{K) is the Fourier transform of the space-time dependent non-abelian field 
A°->^[x) with K = ki + ^2- T^^j are the generators in the fundamental representation, i,j are 
the color indices for the quarks and a is the color index of the non-abelian field. 
Following the above procedure and choosing SU(3) gauge group one obtains 

= tr[T'^T']^ j^^^^^^ d'K (1 - (1 + [K . A%K)K ■ A\-K) 

-K^A\K)- A\-K)]. (9) 
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In the above expression = Using tr\T°'T^] = ^5"''' one obtains 

< = T L>.,, (1 - ^)"V + ^) [A' . A'(K)K . A'(-K) 

-K^A^'iK) ■ A^'i-K)]. (10) 

This is the probabihty for the production of a qq pair from a space-time dependent chro- 
mofield to the lowest order in the couphng constant a^. It can be noted that Eq. (0) is 
valid only for a single flavor of quarks. This equation is simillar to Eq. (^ in the abelian 
theory. 



III. PROBABILITY FOR THE PROCESS A ^ GG 

Now we proceed to compute the probability for the process Ad — * gg via vacuum polar- 
ization. The amplitude for this process (see Fig. 1(a)) is given by 

M = e'^{k,)V^^^liK, h, h)A'^^{K)e'\h). (11) 

In the above expression fci, ^2 are the four momenta of the produced gluons, A°'^{K) is the 
Fourier transform of the non-abelian field A°'^{x) and V^^'{{K^ ki, ^2) is the vertex involving 
a single classical field and two gluons with K = ki + In the conventional method of 
QCD such vertices which involve classical field and gluons are not calculated. Hence, the 
gluon production from a space-time dependent chromofield is not available. However, such 
a calculation is possible in the background field method of QCD which was introduced by 
DeWitt and 't Hooft |jl2|,|l3|. This is because the Feynman diagrams involving a classical 



chromofield and gluons are obtained in the background field method of QCD. We use the 
Feynman rules obtained by the background field method of QCD to compute the probability 
for the process A^ — >■ gg via vacuum polarization. First of all we will briefly describe the 
main differences between conventional QCD and the background field method of QCD before 
studying the above process. 

In the conventional theory of QCD the generating function is 

Z[J] = f [dA] detMc exp{i[S[A] - -G + J-A]), (12) 

where we have not included the quark part for simplicity. In the above expression the gauge 
field action is 

S[A] = -] fd^x {F^f^'f (13) 
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with 

pafj,u _ QfJ-j^au _ Qi^ j^a^i _|_ ^ yafec j^bfi ^cu ^-^^-^ 

for the group with structure constants f"''"^. The other two terms are 

J.A = J d*x r^A''^' (15) 
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and 

G-G = J (fxG^G" (16) 
with G"^ being the gauge fixing term. A typical choice for the gauge fixing term in QCD is 

G" = df.A''^'. (17) 
The matrix element of {Mg{x, y)) is given by 

(*fe(-.^)r' = ^ (18) 

where ^^se^ly)^ derivative of the gauge fixing term under the infinitesimal gauge trans- 

formation 

SA"'' = -P^" e^A^^ + - d^'O". (19) 
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The detMc in the generating functional is written as functional integral over the Faddeev- 
Popov ghost field We mention here that any physical quantity calculated in this method 
is gauge invariant and independent of the particular gauge chosen. However, there is still a 
problem of gauge invariance with some of the quantities like off-shell Green's functions or 
divergent counter terms. This problem arises because in order to quantize the theory one 
must fix a gauge. This means that the total lagrangian we actually use in conventional QCD, 
consisting of the classical lagrangian (which is explicitly gauge invariant) plus a gauge-fixing 
and ghost terms, is not gauge invariant. The background field method is a technique which 
allows one to fix a gauge in quantizing the theory without losing explicit gauge invariance 
which is present at the classical level of the gauge field theory. In the background field ap- 
proach, one arranges things so that explicit gauge invariance is still present once gauge fixing 
and ghost terms are added. In this way Green's functions obey the naive Ward identities of 
gauge invariance and the unphysical quantities like divergent counterterms becomes gauge 
invariant. However, for our purpose of gluon productions from the space-time dependent 
chromofield to the leading order we do not need to address all these issues here. The details 
about the explicit gauge invariance of QCD in the presence of background chromofield can be 
found in []T^ -[15| . In our study we require Feynman rules involving the classical chromofield 
and gluons. For this reason we briefiy outline the procedure of background field method in 
QCD and present the generating functional which generates the Feynman rules involving 
gluons, ghosts and the classical chromofields. 

In the background field method of QCD the gauge fixing term is given by [1^ 



Qa ^ Qt^^ai^ ^ g jabc A^I^Af, (20) 

which depends on the classical background field The variable of the integration in the 



functional integral is the quantum gauge field A"^^ and, following 't Hooft ||T3| the background 
field is not coupled to the external source J . The generating functional depends on J and 
Ad and is given by 
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Z[J, Aa] = J [dAg] detMc exp{i[S[Ag + A,i] - ■ G + J ■ A,]). (21) 
The matrix element of Mq is given by 

(MG(x,y)) =^^^ (22) 

where ^'^^b'^ is the derivative of the gauge fixing term under the infinitesimal gauge trans- 
formation 

^Al^ = -r^^ 9\Al^' + + ^ a^r. (23) 

Writing detMc as functional integral over the ghost field one finds 

Z[J, A,i,^, r] = / [dA,] [dx] [dx*] exp{i[S[Ag + Aa] + S^Host - ■ G + J ■ A, 

+X*^ + Cx])- (24) 

where ^ and ^* are source functions for the ghosts and 

S,host = -Jd'xxiP'6^'' - g'd ,r'\Al^ + A^^) + gt'^A^^d, 

+ 9'nr''Al,^{At^ + Af)]x,. (25) 



The Feynman rules for QCD in the presence of a classical background chromofield 
are constructed from the generating functional which is given by Eq. (^) with the gauge 
fixing term given by Eq. (pO|). We use the Feynman rules obtained in the background field 
method to compute the probability for the process A^ gg. Within this method the vertex 
V^I^^{K, ki,k2), involving a single classical chromofield and two gluons which appear in the 
amplitude (See Fig. 1(a) and Eq. ([TT|)) in Feynman gauge is given by 

V^!;i{K, h, h) = gr'" [2g^.Kx - 2g,xK, - g^xk,], (26) 

where K = ki + k2 and k = ki — k2- The above vertex is different from the three gluon 
vertex usually used in conventional QCD. 
Using Eq. (0) in (|) we find 

Tgi = h:,,,^A^^{K)A^'^'{-K)V^^^^,{K, k)VX!''^\K, k)e'''{k,)e*''''' {k,)e^\k2)e*^''' (k^), (27) 

where we have used A*{K) = A{—K). The factor | is the weight factor given in the Feynman 
amplitude in order to obtain a correct gauge field and gauge-fixing action. To obtain correct 
and physical results we have to use the appropriate projection operators for the transverse 
polarization states of the gluons. For this purpose we proceed as follows. First of all we use 
the polarization sum 

and then substract the corresponding ghost contributions (see Fig. 1(b)). 
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Using Eq. (^), (^) and performing the calculation we find 

Tgi = -^A^^^{K)A-'>^\-K)V^^l{K, k)VX^^,{K, k)g^^'g'^' = Ng' [Ak, ■ k^A\K) ■ A\-K) 

-3ki ■ A\K)k2 ■ A''{-K) - 3k2 ■ A''{K)ki ■ A%-K) - ki ■ A%K)ki ■ 
-k2-A%K)k2-A''{-K)], (29) 

where we have used f°-^<^ f°''''<^ = N6'^^' in SU(N) gauge group. 

Now we proceed to compute the probability for the process Ad — ^ gg. To obtain the 
probability for the above process we require that = {ki + fcs) > with K'^ > 0. We 
recall that for real gluons A;^ = A;| = 0. With the above requirements we proceed to perform 
the integration in Eq. (|^) using Tgi from Eq. (^9]). As gluons are similar particles we 
multiply a factor | in the phase space and find 

^(1) 1 Ng^ f j4r^/i/r^0N fd^kid?k2 



gi 



4 327r2 7x2>o 



[Sg^^'K'^ - SKf^K" + Ak^'k"] (30) 



where k = ki — k2. 

Now, we evaluate the corresponding ghost diagrams which have to be substracted from 
the above result. The amplitude for the ghost part (see Fig. If,) is given by 



M'g'^ = g^k^Al. (31) 



so that 



T,, = \mI^^ mX (32) 



Using Eq. (R^) in (Rl) we find 



"'S' = iiiS U/" '"^'■°> / ^^*'"<^^' - - h)A:(i<)A:{-i<)ikV] (33) 

Substracting Eq. (|33D from (|30|) we find 

[8c/^^ - SRf'K'' + 2k^'k''] (34) 
To perform the integration in the above equation we proceed as follows. First of all using 

d'k2 



d^k2 5{kl) e{kl) (35) 



we check that 



/ ^^^^^^'Hi^ -h- k2)g^^ = 2ng^r (36) 

J /Ci /Co 
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Using Eq. (|36D we find 



.0 .0 -k,- k,)k^,k^, = f ^^6^'\K -k,- k,)ki^k\ 



-[i^V + 2X''^1, (37) 
6 



and 



d?kid?'k2 ^(A\, , , rd^kid^k2 



1 2 1 2 



5(4)(x - A;i - A;2)A;5'A;r = / ^^S^'\K -h- k2)k!^kl^ 



'^[-K^gt^^ + AK^'K'']. (38) 



With the help of the above relations the integration in Eq. (|^) can be easily performed. 
Using (|36|), (|37|) and ( pSD we obtain from Eq. (0) 

^ (fK e{K^) [K^A''{K) ■ A%-K) - K ■ A%K)K ■ A^^-K)] (39) 



Agg 24 

In the above expression the repeated indices are summed from a = 1, ...{N"^ — 1). For SU(3) 
gauge group we obtain 



W 



(1) 



Agg 



I d^K e{K^) \K^A''{K) ■ A^{-K) - K ■ A''{K)K ■ A\-K)] (40) 



This is the probability computed from the process A ^ gg via vacuum polarization. It can 
be seen that the above expression is transverse with respect to the momentum of the field 
K, i.e, when A'^ is replaced by we get W^gg = 0. This result is a part of the expression 
for the total probability (including higher order terms in gA) for the production of gluon 
pairs. 



IV. CONCLUSIONS 

Using the background field method of QCD, we have computed the probability for the 
processes A qq, gg via vacuum polarization, in the presence of a space-time dependent 
chromofield A. These processes are similar to A e^e~ in QED. For any arbitary classical 
non-abelian chromofield one has to check the gauge invariance of the result with respect to 
the general non-abelian local gauge transformation which is known as type I gauge trans- 



formation [|T6[. Under this gauge transformation the classical chromofield transform like: 

A!:,^UA^^,U-'-'-{d,U)U-\ (41) 

The gluonic and ghost fields transform like: 

A^^UA>'qU-\ (42) 

and 
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X' ^ UxU-^ (43) 

respectively. 

In these situations one has to add the results of the higher order processes (for the 
gluon pair production) to the leading order results obtained in this paper. The gluon pair 
production amplitude which is obtained from the first order S matrix contains higher order 
terms in gA. This is due to the non-abelian nature of the gluonic and classical chromofield. 
The result of the gluon pair production probability which contains both leading order and 



higher order terms in gA (but first order in S matrix) will be reported elsewhere |jT^. In this 
paper we have presented only the result from the leading order processes A — > gg. This result 
is a part of the expression for the total probability of the gluon production from a space- 
time dependent chromofield [|l7l. Quark and gluon production from a space-time dependent 
chromofield will play a crucial role in the production and equilibration of the quark-gluon 
plasma in ultra relativistic heavy-ion collsions. Our main goal is to study the production 
and equilibration of the quark-gluon plasma by solving relativistic non-abelian transport 
equations for partons (see Eq. ([^)) with parton production from a space-time dependent 
chromofield taken into account. Such work requires extensive numerical computation |T 
and needs a separate publication. 
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Fig. 1 Lowest order vacuum polarization diagrams. Fig. 1(a) is the 
gluon diagram and Fig. 1(b) is the corresponding ghost diagram. 
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